Effect of interface bonding on spin-dependent tunneling from the oxidized Co surface 
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We demonstrate that the factorization of the tunneling transmission into the product of two sur- 
face transmission functions and a vacuum decay factor allows one to generalize Julliere's formula 
and explain the meaning of the "tunneling density of states" in some limiting cases. Using this 
factorization we calculate spin-dependent tunneling from clean and oxidized fee Co surfaces through 
vacuum into Al using the principal-layer Green's function approach. We demonstrate that a mono- 
layer of oxygen on the Co (111) surface creates a spin-filter effect due to the Co-0 bonding which 
produces an additional tunneling barrier in the minority-spin channel. This changes the minority- 
spin dominated conductance for the clean Co surface into a majority spin dominated conductance 
for the oxidized Co surface. 



I. INTRODUCTION 

Spin-dcpcndcnt tunneling (SDT) in magnetic tunnel 
junctions (MTJs) is a dynamically developing area of re- 
search that attracted a lot of attention due to promising 
applications in non- volatile random access memories and 
next-generation magnetic field sensors (for a recent re- 
view of SDT sec Ref.Q). The experimental efforts have 
succeeded in achieving large reproducible tunneling mag- 
netoresistance (TMR) in MTJs^ but also raised funda- 
mental questions regarding the nature of SDT. One such 
question is the role of the ferromagnet/insulator inter- 
faces in controlling the spin polarization (SP) of the tun- 
neling conductance. The SP is usually defined as P = 
(G T - G L )/(G 1 + G x ), where G a = (e 2 /h) £ k|| T CT (k,|) is 
the conductance for spin channel a, T a is the transmis- 
sion function, and k|| is the lateral wave vector. 

Commonly the expected spin dependence of the tun- 
neling current is deduced by considering the symmetry 
of the Bloch states in the bulk ferromagnetic electrodes 
and the complex band structure of the insulatori^4 By 
identifying those bands in the electrodes that are allowed 
by symmetry to couple efficiently to the evanescent states 
decaying most slowly in the barrier one can predict the 
SP of the conductance. However, this approach has two 
deficiencies. First, it assumes that the barrier is suffi- 
ciently thick so that only a small focused region of the 
surface Brillouin zone (SBZ) contributes to the tunneling 
current. For realistic MTJs with a barrier thickness of 
about 1 nm this assumption is usually unjustified. Sec- 
ond, symmetry considerations alone applied to bulk ma- 
terials are not always sufficient to predict the SP. It is 
critical to take into account the electronic structure of 
the ferromagnet/barrier interfaces which, as it was shown 
both experimentally^ and theoretically^ controls SDT. 

An important mechanism by which the interfaces affect 
the SP of the tunneling current is the bonding between 
the ferromagnetic electrodes and the insulator^ This ef- 
fect was put forward to explain positive and negative val- 



ues of TMR depending on the applied voltage in MTJs 
with Ta2C>5 and T^Os/A^Os barriers^ and to elucidate 
the inversion of the SP in Co/SrTi03/Lao.67Sro.33MnC>3 
MTJsiS So far there are no theoretical studies explaining 
the microscopic origin of this phenomenon. 

In this paper we report the results of first-principles 
calculations of SDT from clean and oxidized Co surfaces 
through vacuum into Al and demonstrate the crucial role 
of the bonding between Co and O atoms. This system 
was chosen for investigation because the MTJs based 
on alumina have predominantly O-terminated C0/AI2O3 
interfaces^ By replacing alumina by vacuum we can ig- 
nore the complexity of the atomic structure of the amor- 
phous alumina and focus on effects of surface oxidation. 
Moreover, this system can be directly studied using spin- 
polarized STM. [n] 

We show that a monolayer of oxygen on the Co sur- 
face creates a spin-filter effect due to the Co-0 bond- 
ing by producing an additional tunneling barrier in the 
minority-spin channel. This reverses the sign of the SP 
from negative for the clean Co surface to positive for the 
oxidized Co surface, thus revealing the crucial role of in- 
terface bonding in SDT. 



II. SURFACE TRANSMISSION FUNCTIONS 

We approach the tunneling problem in the spirit of 
perturbation theory^ We consider the system consist- 
ing of "left" and "right" leads separated by a relatively 
thick barrier and assume that two-dimensional transla- 
tional periodicity in transverse directions is preserved for 
both electrodes including their surfaces (although it may 
be different for each subsystem) . Every Bloch wave with 
a lateral quasi-wave vector kj^ coming from the left lead 
has a decay tail in the vacuum composed of the waves 
with transverse wave vectors kj^ + G, where G^ are the 

reciprocal lattice vectors of the SBZ of the left leadii 3 - 
At sufficient distances from the surface (typically just a 
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few lattice parameters for low-index surfaces) all waves 
with Gj ^ vanish and may be neglected in the Fourier 
expansion of the Bloch wave in the vacuum. On the 
other hand, the quasi-wave vector kj^ is the good quan- 
tum number for the right lead. This means that k|| is con- 
served across the entire system even if there is no common 
in-plane periodicity and k|% kj^ are defined with respect 
to different SBZs. In this case, each tunneling eigenstate 
is almost identical to an evanescent plane wave in the 
central region of the barrier. 

For a given k|| the transmission function is the sum of 
the transmission coefficients for all tunneling eigenstates 
corresponding to all incoming Bloch waves with this k|| 
in the left electrode. At the same time, each transmission 
coefficient for a given eigenstate contains a sum over out- 
going states in the right electrode with the same k|| . Let 
us choose a reference plane in the vacuum region at a suf- 
ficient distance from the surface of an electrode, so that 
the eigenstates for all kj| are already indistinguishable 
from the barrier eigenstates at this plane (see Fig. QJ. 
For each tunneling eigenstate the amplitude of the bar- 
rier eigenstate between the reference planes is the only 
parameter coupling the left and right electrodes. Then, 
the S-matrix element coupling the states in the two elec- 
trodes may be written as 
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FIG. 1: Geometry of a tunnel junction with a vacuum barrier. 
The top graph schematically shows the potential barrier for 
electrons at the Fermi level. Dashed lines show the positions 
of the reference planes assigned to the electrodes for the calcu- 
lation of the surface transmission functions. Each scattering 
state is described by four amplitudes: A v of the incoming 
Bloch wave; B r , B' r of the surviving evanescent wave at the 
reference planes assigned to the left and right electrodes; and 
C q of the transmitted Bloch wave in the right electrode (there 
may be several transmitted waves). The bottom scheme de- 
picts atomic layers in the electrodes and their labelling used 
in the text. 
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where p is the incoming and q the outgoing Bloch states 
in the left and right leads respectively, r and r' denote 
the same vacuum eigenstate at the left and right refer- 
ence planes. We omitted the dependence on k|| for all the 
S-matrices in Eq. for brevity. The vacuum S-matrix 
S rr i simply describes the exponential decay of the wave 
function in the vacuum. Note that no summation is im- 
plied in |JI)|. because the state r is uniquely defined by 
k|| . The simple product of S-matrices in Eq. without 
any multiple scattering terms is a consequence of our as- 
sumption that the barrier is sufficiently thick. Thus, we 
see that the transmission function r(kii) of the MTJ is 
factorized: 



T a (k|,)=*J(k||)exp[-2K(k||)d]tS(k||). 



(2) 



Here we replaced |<SV r '| 2 by its explicit exponential form 
with 
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where <f> is the work function, and d is the distance be- 
tween the reference planes assigned to the electrodes as 
shown in Fig. ^ All the information about the prop- 
erties of individual surfaces is described by the surface 
transmission functions (STF) tf 



fcO<ii)=E x 1 > **( k n) = £ 



C a 
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where the four amplitudes characterize the behavior of 
tunneling eigenstates at the two surfaces (see Fig. QJ. 
The different definitions of and tn are due to the fact 
that they obey different boundary conditions. Specif- 
ically, tjj and tn are obtained by matching the Bloch 
wave functions with the vacuum eigenstates, decaying 
and growing into the vacuum respectively. The defini- 
tion of ti, implies the solution of a scattering problem 
for the incoming wave with amplitude A p , and B r is the 
amplitude of this scattering eigenstate at the reference 
plane in the vacuum. On the other hand, tn describes 
an inverse scattering problem in which the exponentially 
decaying wave in the vacuum with amplitude B r i at the 
reference plane is scattered on the right surface; here C q 
is the amplitude of the outgoing Bloch state q in the right 
electrode for this "eigenstate." Note that physically, this 
state is forbidden because it grows to infinity in the vac- 
uum, but it is still a formal solution of the Schrodinger 
equation with the specified boundary condition at infin- 
ity. 

Each surface may be considered as forming the left or 
the right interface of the tunnel junction. It is straight- 
forward to show that, just as in the case of a transmis- 
sion matrix connecting the propagating states in the two 
electrodes^ the requirement of flux conservation results 
in the reciprocity condition tz,(kii) = ^(kii) = i(kii) 
for any (laterally periodic) surface, as long as the ap- 
propriate normalization of the wave functions is chosen. 
Specifically, all Bloch waves in the electrodes should be 
normalized for unit flux, while the vacuum eigenstates 
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should be normalized for unit "imaginary flux" n/m. 
The easiest way to establish this reciprocity condition 
is to prove it for free electrons in a semi-infinite poten- 
tial well, and then to use this free-electron system as a 
second electrode in a junction. Since the total transmis- 
sion function for a two-electrode system is reciprocal, 14 
this proves the above reciprocity condition for the STF. 
Note that this proof does not rely on time reversal sym- 
metry because time reversal also replaces kn by — kii. 
However, this symmetry gives another useful relation, 
i(k|| , H) = £(— kii , H) where H is the external magnetic 
field. 

It is easy to see from Eq. Q with the chosen unit flux 
normalization that the STF is proportional to the Fermi- 
level value of the k||- and energy-resolved charge density 
for the given spin, which is generated by the incoming 
Bloch states and taken at the reference plane (any lo- 
calized surface states are excluded). Indeed, the energy- 
resolved charge density may be written as 



where the k||- and energy-resolved charge density is 



r (r,k||,£)=^Lf k||P (r; 



(5) 



(6) 



where the summation is over the incoming Bloch states 
with the given k|| and _E, and the superscript F empha- 
sizes that the scattering eigenstates i>^ llP are normalized 
so that the incoming Bloch waves carry unit flux normal 
to the interface. Setting A p = 1 in Eq. (@J, we see that 
t(k|| ) is given by Eq. © where r is taken at the reference 
plane (since we neglect all components with Gi ^ at 
the reference plane, p(r,\s.«,E) does not depend on r||, 
the location within this plane). 

Although we considered a vacuum tunneling barrier, 
the analysis can be extended to other physically impor- 
tant cases of insulating barriers. Indeed, the main re- 
quirement for the validity of Eqs. ill 121) is that the com- 
plex band structure be predominantly represented by 
a single evanescent wave for each k|| . At a minimum, 
this premise must hold for the "active" regions of the 
SBZ that contribute appreciably to the tunneling cur- 
rent. This means that other tunneling states must have 
a notably larger imaginary part of the wave vector com- 
pared to the dominant one. In the case of a vacuum 
barrier, this amounts to the neglect of all vacuum eigen- 
states with Gi 7^ 0. In practice this criterion is well sat- 
isfied for sufficiently thick wide-gap tunnel barriers, such 
as sp-bonded oxides. However, for any particular barrier, 
this assumption has to be carefully verified by checking 
the complex band structure of the barrier for the pres- 
ence of additional slowly decaying states in the active 
regions of the SBZ. It is important to note that the area 
of the active region quickly shrinks as the barrier thick- 
ness is increased^ indicating that the single evanescent 



state criterion will often lead to an additional require- 
ment of "sufficient thickness" of the barrier. Note that 
this requirement is not related to a similar one stemming 
from the neglect of back-scattering. 

The factorization J5J allows us to study tunneling be- 
tween completely different electrodes. For practical pur- 
poses, the STF can be calculated using T CT (ki|) obtained 
for a symmetric MTJ by factoring out the vacuum decay 
factor for the given choice of reference planes, and then 
taking the square root. The resulting STFs for different 
surfaces may then be convolved with an appropriate vac- 
uum decay factor, exp(— 2nd), to obtain the transmission 
functions for asymmetric MTJs. In particular, T^kn) for 
the antiparallel (AP) magnetization of the two electrodes 
(and the TMR) can be found from the up- and down-spin 
STFs calculated from the transmission function for the 
parallel (P) magnetization. 



III. GENERALIZATION OF JULLIERE'S 
FORMULA 

Tunneling magnetoresistance is often discussed in 
terms of Julliere's formula 15 



TMR: 



2P l Pr 
- - PlPr 



(7) 



where Pl , Pr are the "spin polarizations" of the left and 
right electrodes. Eq. Q is derived from the empirical 
expression G a oc p a L p R , where p a Ll p R are the "tunneling 
densities of states" of the electrodes, and the spin po- 
larization is defined as P = (p^ — p l )/(/0 T + p^)- The 
popularity of this formula is due to the fact that it usu- 
ally agrees reasonably well with experiment if the spin 
polarizations, which are directly related to the "tunnel- 
ing densities of states," are taken from Meservey-Tedrow 
experiments 16 with the same barrier. However, the valid- 
ity of Julliere's formula has been debated for a long time, 
and the reasons for its apparent agreement with experi- 
ment are unclear. The physical meaning of the "tunnel- 
ing density of states" is also unclear, but it is obvious 
both from elementary quantum mechanics and from ex- 
periments that the tunneling properties of a magnetic 
heterostructure are determined not by the ferromagnet 
alone, but rather by the ferromagnet/barrier combina- 
tion and by the structure of the interface. A number 
of explicit first-principles calculations for idealized MTJs 
without disorder (see, e.g., Ref. \\% confirmed this fact. 
However, it was suggested 1 ^ that phase decoherence due 
to disorder which is always present in realistic MTJs may 
recover the factorization of the tunneling conductance in 
a product of "transport densities of states," which are 
essentially equal to the regular densities of states at the 
surfaces of the electrodes if there are no resonant localized 
states in the barrier. Moreover, it was shown 1 ^ within a 
single-band tight binding model that in the limit of strong 
disorder one recovers Julliere's formula J7J) by identify- 
ing Pl , Pr with the measurable spin polarizations of the 
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tunneling current for the same electrode/barrier systems. 
Therefore, it seems that there are good reasons for the 
widespread use of Julliere's formula, and it is highly de- 
sirable to elucidate these reasons. 

Let us explore the connection between Eq. and 
Julliere's formula JJJ. In Eq. (J2J) the simple product of 
the "tunneling densities of states" is replaced by a con- 
volution of STFs, which explicitly include the effects of 
bulk densities of states and of the surface structure. Thus 
we can consider Eq. J5J) as a generalization of Julliere's 
formula for an ideal MTJ with no disorder. 

However, we may go further and identify limiting cases 
where Eq. (0 can be directly related to Julliere's for- 
mula, providing formal definitions of the "tunneling den- 
sities of states" of the electrodes appropriate for these 
cases. First, consider the case of a disordered insulat- 
ing barrier. Such a barrier may be characterized by its 
eigenstates, half of which are decaying from left to right, 
and the other half from right to left. Although these 
eigenstates do not have a conserved kn any more, it is 
still clear that tunneling will be dominated by Feynman 
paths that do not "loop back," because each path carries 
a weight decaying exponentially with its length (see be- 
low). Therefore, we may still write an expression similar 
to Eq. (|TJ neglecting back-scattering, but now we should 
sum up over all barrier eigenstates (now defined in real 
space). Within this formulation S rr ' is still diagonal be- 
cause it describes the decay of a single eigenstate. 

The weight of a Feynman path in the imaginary-time 
functional integral often used for tunneling problems 20 
(with Euclidean action written in its reduced Mauper- 
tuis form; see, e.g., Ref. l2l|) is given by exp [— J n(l)dl] 
up to a prefactor, where the integral is taken along the 
path, k = [2m(V — E)] 1 / 2 and V(r) is the potential. In an 
ordered insulator many paths with similar weights con- 
tribute to the path integral resulting in the formation of 
the complex band structure. However, in a disordered 
insulator the tunneling current may be dominated by 
Feynman paths running close to a relatively small num- 
ber of "easy" paths with locally maximum weights, i.e. 
by imaginary-time classical paths ,22*21 If there is only 
one such channel or one class of channels with similar 
properties (e.g., due to surface roughness), Eq. J3J) will 
produce Julliere's formula where p a is simply the Fermi- 
level value of the energy-resolved charge density given by 
Eq. @ integrated over k|| and taken at some reference 
point within the channel. (Now each term describes the 
scattering eigenstate corresponding to the single incom- 
ing Bloch wave with the given k|| .) Like STF, this quan- 
tity does not depend on the properties of the other side 
of the barrier. This conclusion agrees with the results 
of Ref. showing that the tunneling current through a 
strongly disordered barrier is dominated by a small num- 
ber of random configurations, and that Julliere's formula 
is also recovered in this limit. 

Now consider the case when disorder is weak close to 
the interfaces, but remains strong in the insulator. Ob- 
viously, the S-matrix of the disordered insulator in kn 



representation will be essentially a random matrix, and 
after averaging Eq. @ thus yields Julliere's formula with 
p a oc Ylk ^(^il)- This case is the easiest from the com- 
putational point of view, because the STFs may be di- 
rectly calculated for a k|| -conserving MTJ. 

It is instructive to compare this result with the con- 
clusions of Mathon and Umersk:- 1 — on the applicability of 
Julliere's formula obtained using the transfer Hamilto- 
nian formalism. Our approach shares in common with 
Ref. the neglect of multiple reflections across the 
junction. However, the assumption of constant matrix 
elements (hopping integrals) for all Bloch waves made 
in Ref. [Tj| completely removes all physical effects con- 
nected with orbital- and spin-dependent bonding at inter- 
faces. This obviously contradicts the experimental find- 
ings showing that the spin polarization of the tunneling 
current and magnetoresistance strongly depend on the 
type of barrier usedi In our approach, the STFs for the 
electrodes allow us to encapsulate the effects of the inter- 
face structure and provide the proper dependence of the 
tunneling current on barrier type. Julliere's formula ob- 
tained in the limiting case of full dccohcrence inside the 
insulator is expressed in terms of the spin polarization 
actually measured in the Meservey-Tedrow experiment 
(assuming that the superconductor acts as an ideal, non- 
biased spin detector). 

Finally, for very thick k|| -conserving barriers the tun- 
neling current may be carried predominantly by a close 
vicinity of some special k|| -points in the SBZ (e.g., the 
r point). In this case, the tunneling density of states is 
simply equal to the value of T(k||) at this k||. 

It is not clear a priori whether any one of these three 
limiting cases is directly applicable to realistic MTJs, al- 
though it seems that disorder in the insulator together 
with the "channelization" of the tunneling current are 
both likely to play a major role. However, the emergence 
of Julliere's formula in these different scenarios suggests 
that it may actually have a rather wide range of appli- 
cability. In general, the "tunneling density of states" 
should be identified with some appropriately averaged 
energy-resolved charge density taken at the Fermi level 
at a sufficient distance from the interface within the bar- 
rier. Unlike the bulk density of states, this function fully 
takes into account the relevant properties of the surface. 



IV. TUNNELING FROM CLEAN AND 
OXIDIZED Co (111) SURFACES THROUGH 
VACUUM INTO Al 

We calculated the transmission functions using the 
principal-layer Green's function approach^ based on 
the tight-binding linear muffin-tin orbital method (TB- 
LMTO) in the atomic sphere approximation (ASA) and 
the transmission matrix formulation of Ref. |2j All 
the atomic potentials were determined self-consistently 
within the supercell approach using the TB-LMTO-ASA 
method. The vacuum barrier was modeled using empty 



spheres in the positions corresponding to the continu- 
ation of the crystal lattice of the electrodes. We have 
also performed full-potential LMTO calculations 24 which 
confirmed all main features of the band structure of the 
oxidized Co (111) surface discussed below. 

We checked the validity of factorization by calculat- 
ing To^kn) for (100) and (lll)-oriented fee Co electrodes 
with parallel magnetizations, taking the square root, and 
convolving tf (kn ) with ij_ (kii ) . Then, the result was com- 
pared with the independent calculation for the antipar- 
allel configuration in a range of energies. The agreement 
was always excellent (better than 1%), except for a cou- 
ple of specific energies for a (100) MTJ with four vacuum 
"monolayers" where narrow resonances appear in the mi- 
nority channel^ If the vacuum barrier is extended to 8 
monolayers, excellent agreement is restored. 

Using the factorization J2J) we investigated the SP of 
the conductance from ferromagnetic electrodes to a non- 
magnetic material, Al (111), which served as a detector 
of the tunneling SP in the spirit of the Meservey-Tedrow 
experiments^ As expected, the calculated STF of Al is 
free-electron-like, having almost perfect Gaussian shape 
originating from the vacuum decay factor up to the ref- 
erence plane. Therefore, this surface may be considered 
as equally transparent for all Bloch waves, and the total 
transmission function for a MTJ with Al spin-detector 
electrode is essentially a product of the other electrode's 
STF and the vacuum decay factor. 

First, we discuss the properties of a Co/vacuum/ Al 
MTJ with a clean Co (111) surface. Figs. [2^, b show the 
k|| -resolved transmission for the majority- and minority- 
spin electrons within the SBZ of Co (111). The Fermi 
surface of Co viewed along the [111] direction has holes 
close to the T point with no bulk states in both spin 
channels, which results in zero conductance in this area. 
The majority-spin transmission (Fig. varies rather 
smoothly and is appreciable over a relatively large area 
of the SBZ. On the other hand, the minority-spin trans- 
mission (Fig. ^jp) has a narrow crown-shaped "hot ring" 
around the edge of the Fermi surface hole. The analysis 
of layer and k|| -resolved density of states (DOS) shows 
that it is not associated with surface states^ but rather 
with an enhancement of bulk k|| -resolved DOS near the 
Fermi surface edge (compare Fig. [2) with Fig. 0J;). 

As seen from Figs. |2Kjb, the Fermi surface hole is 
smaller for majority spins. Therefore, the conductance 
should become fully majority-spin polarized in the limit 
of very thick barriers. However, since the Fermi surface 
hole is also quite narrow for minority spins, positive SP is 
only achieved at very large barrier thicknesses d ~ 10 nm, 
while for typical values of d ~ 2 nm the SP is about —60% 
and depends weakly on d. 

The oxidized Co surface was modeled by an O mono- 
layer placed on top of the Co(lll) electrode. The equilib- 
rium atomic structure of this surface was found using the 
pseudopotential plane-wave method^ within the gener- 
alized gradient approximation. We used both types of 
stacking: ABC A and ABCB, where the last symbol des- 
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FIG. 2: k||-resolved transmission (logarithmic scale) from 
clean and oxidized (111) Co surfaces through vacuum into 
Al. (a) Clean surface, majority spin, (b) Clean surface, mi- 
nority spin, (c) Oxidized surface, majority spin, (d) Oxidized 
surface, minority spin. The vacuum layer thickness is 2 nm 
for clean and 1.7 nm for oxidized Co surface. The first sur- 
face Brillouin zone is shown. Units are 10 -11 for (a), (b) and 
1(T 14 for (c), (d). 



ignates the position of the O monolayer. The O atoms are 
assumed to lie in symmetric positions above the second 
(S2) or first (SI) sub-surface Co layer, respectively (these 
layers are shown in Fig. ^in the absence of the O mono- 
layer). We allowed the O layer and two Co layers (SI and 
S2) to relax in the direction normal to the surface, while 
the positions of atoms in deeper layers (S3,. . . ) were kept 
fixed. The energies of both equilibrium configurations of 
oxygen monolayers were found to be very close to each 
other. All results of interest in the present context are 
quite similar for these two stackings. Below all specific 
data are given for the ABCB stacking. The equilibrium 
interlayer distances were found to be 2.14 A between lay- 
ers S3 and S2, 2.18 A between S2 and SI, and 1.08 A be- 
tween SI and O layers, compared to 2.07 A between the 
adjacent Co layers in the bulk. The Co-0 bond length is 
equal to 1.82 A. 

Presence of oxygen at the surface of cobalt raises the 
question of whether electron correlations similar to those 
characteristic for transition-metal oxides may be strong 
enough to induce significant changes in the band struc- 
ture at the oxidized Co surface. However, the enhance- 
ment of correlations in oxides is due to much weaker 
screening of Coulomb interaction compared to the metal- 
lic state. On the other hand, cobalt atoms below the oxy- 
gen monolayer preserve the close-packed configuration of 
bulk cobalt except for the top three nearest neighbors 
being absent out of twelve. Therefore, it is reasonable to 
expect that screening of Coulomb interaction in the 3d 
shell is not much weaker compared to the bulk. For this 
reason, we believe that LDA electronic structure of the 
oxidized Co surface is correct as far as the main features 
are concerned. 
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FIG. 3: Band dispersions along the TM direction for (a) ma- 
jority and (b) minority spin electrons. Energy is referenced 
from Ef- The bonding Co-0 surface bands are marked B, 
and the pure antibonding surface band is marked A. 



The oxygen monolayer dramatically changes the elec- 
tronic structure of the underlying Co layer making this 
layer almost magnetically-dead. This change can be un- 
derstood from band dispersion plots shown in Fig. [21 For 
each spin, the free-standing oxygen monolayer has three 
energy bands deriving from 2p states, each doubly degen- 
erate due to a z reflection symmetry (z is the axis normal 
to the surface). When the monolayer is deposited onto 
the Co surface, the degeneracy is lifted, and two sets of 
three bands each are formed corresponding to bonding 
and antibonding mixing of oxygen and cobalt orbitals. 
The three bonding bands (marked B in Fig.[3J) lie well be- 
low the bulk Co 3d band, whereas the antibonding states 
are close to the Fermi energy Ep. As a result of this 
bonding the local DOS for the SI layer at Ep is strongly 
reduced, so that, according to the Stoner criterion, mag- 
netism in this layer is almost completely suppressed. The 
magnetic moment of Co atoms in the SI layer is only 0.17 

Transmission of propagating bulk states from the elec- 
trode through the barrier is very sensitive to the degree 
of mixing of these states with the antibonding surface 
states. This mixing is controlled by a selection rule 
which follows from the fact that all bands can be clas- 
sified as "even" or "odd" according to their symmetry 
with respect to TM reflection (TM is the projection of 
a reflection plane). Although this classification is exact 
only along the TM direction, it is approximately valid 
throughout the entire SBZ. 

According to this classification, two of the three surface 
bands are even, and one is odd. On the other hand, the 
free-electron-like band of bulk Co which forms the only 
Fermi surface sheet for majority-spin electrons is even, 
while the minority-spin states on the Fermi surface sheet 
closest to the f point are odd. Even and odd bands are 
orthogonal and can not mix. 




FIG. 4: ky-resolved local DOS at Ef (arbitrary units) for 
the oxidized Co surface: (a) S6 layer, majority spin; (b) SI 
layer, majority spin; (c) S6 layer, minority spin; (d) SI layer, 
minority spin. 

This selection rule results in the principal difference 
between the majority- and minority-spin transmission. 
The majority-spin bands are shown in Fig.|3Ji. One even 
and one odd antibonding surface bands (marked A) are 
degenerate at the T point (at about 0.8 eV above Ef). 
At a short distance from the T point both bands enter 
the continuum of bulk states. The odd band does not 
mix with the bulk states along the TM line and remains 
almost flat due to repulsion from a lower-lying band. On 
the other hand, the even band readily mixes with the 
free-electron-like majority-spin band crossing the Fermi 
level and completely loses its surface localization. This 
is evident from Figs. 2Ji,b which show the k||-resolved 
DOS of majority-spin electrons for the bulk-like S6 and 
surface SI Co layers. The DOS for the SI layer (Fig. Ha) 
is appreciable in the entire area of the SBZ where bulk 
states are available and has no sharp features that might 
indicate localized surface states. This implies that the 
bulk majority-spin states extend to the very surface of 
the electrode and therefore can readily tunnel through 
the barrier. 

The situation is very different for minority-spin states. 
Although the odd surface band is again almost flat and 
lies above Ep, the even surface band crossing the Fermi 
level does not mix with the odd minority-spin band. As 
a consequence, the k||-resolved DOS for Co(Sl) layer 
(Fig.^Ji) is large only along the curve lying at the periph- 
ery of the SBZ where the Co-0 antibonding surface band 
crosses the Fermi level (oxygen DOS looks very similar). 
As a result, the bulk minority-spin states responsible for 
most tunneling transmission from the clean surface only 
extend up to the S2 layer, encountering a band gap in the 
SI and oxygen layers. Thus, an additional tunneling bar- 
rier is introduced in the minority-spin channel, and the 
total SP of the tunneling transmission becomes almost 
100% positive, which is evident from Figs. |2t,d. 

The predicted effect of interface bonding is not limited 
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only to the Co(lll) surface. We have also calculated 
the transmission from clean and oxidized Co(100) and 
Ni(lll) electrodes and found that surface oxidation re- 
verses the SP due to the bonding between Co or Ni with 
O. As it was shown earlier, the reversal of the SP also 
occurs for the Fe(100) surfaced 

V. CONCLUSION 

We have shown that the problem of calculating of the 
transmission function for a sufficiently thick insulating 
barrier is reduced to the solution of three separate prob- 
lems, namely the penetration of the bulk wave functions 
into the barrier from both sides, and the behavior of 
the evanescent barrier eigenstates. This separation pro- 
vides a natural generalization of Julliere's formula. We 
identified three limiting cases when the original Julliere's 
formula is recovered. The "tunneling density of states" 
in this formula is identified with an appropriately av- 
eraged energy-resolved charge density generated by the 
bulk Bloch states within the barrier and taken at the 
Fermi level. 

Using the factorization of the transmission function 
into a product of surface transmission functions and a 
barrier decay factor we calculated the spin polarization 
of the tunneling current from clean and oxidized Co (111) 
surfaces through vacuum into Al. We showed that the 



bonding between Co and O atoms at the oxidized sur- 
face controls SDT by creating an additional barrier for 
minority-spin electrons, which results in a reversal of the 
spin polarization. 

Experimentally, the reversal of the SP associated with 
surface oxidation may be detected using spin-polarized 
STM measurements, — Since the ferromagnetic tip is sen- 
sitive to the SP of the total local DOS above the sur- 
face (see, e.g., Ref. Esl) . the TMR in the system sur- 
face/vacuum/tip should change sign when the Co surface 
is oxidized. In other words, for the clean Co (111) surface 
the tunneling current should be higher when the magne- 
tizations of the tip and the surface are aligned parallel 
(the dominating minority channel is then open), but for 
the oxidized surface it should be higher for the antipar- 
allel configuration. 
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